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by 

Jillian A. K. Stott 
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Abstract 

We consider the effects of viscosity on the inviscid stability of the Batchelor (1964) vortex 
in a compressible flow. The problem is tackled asymptotically, in the limit of large (stream- 
wise and azimuthal) wavenumbers, together with large Mach numbers. This problem, with 
viscous effects neglected, was discussed in Stott & Duck (1994). The authors found that 
the nature of the solution passes through different regimes as the Mach number increases, 
relative to the wavenumber. This structure persists when viscous effects are included in 
the analysis. In the present study, as in that mentioned above, the mode present in the 
incompressible case ceases to be unstable at high Mach numbers and a centre mode forms, 
whose stability characteristics are determined primarily by conditions close to the vortex 
axis. 

We find generally that viscosity has a stabilising influence on the flow, whilst in the 
case of centre modes, viscous effects become important at much larger Reynolds numbers 
than for the first class of disturbance. 


* This work was partially supported by the National Aeronautics and Space Adminis- 
tration under NASA Contract No. NASl-19480 while the author was in residence at the 
Institute for Computer Applications in Science and Engineering (ICASE), NASA Langley 
Research Center, Hampton, VA, 23681-0001. 
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1. Introduction and formulation 


The work presented in this paper is a follow on study to that of Stott & Duck (1994), where 
the inviscid stability of the Batchelor (1964) trailing-line vortex in a compressible flow was 
considered. The motivation for this research is due to relevance to tornadoes, combustion 
optimisation and also the breakdown of such vortices behind aircraft. 

The first authors to investigate viscous effects in the study of the stability of the 
Batchelor (1964) vortex were Lessen &: Paillet (1974). A numerical approach was adopted 
and it was shown the n = — 1 mode (where n is the azimuthal wavenumber) was the 
most unstable mode. This finding was contrary to previous inviscid studies i.e. Lessen, 
Singh & Paillet (1974), which had indicated that disturbances with — n 1 were the most 
dangerous. However the work of Lessen, Singh &; Paillet (1974) considered only relatively 
low Reynolds numbers (i.e. Re ~ 100), which is unrealistic for comparison with experimental 
flow situations and with inviscid studies. This apparent anomaly remained unresolved until 
the study by Stewartson (1982), in which the problem was approached asymptotically (based 
on the inviscid work of Leibovich & Stewartson 1983) which considered the dual limits of 
large Reynolds numbers and wavenumbers. The results obtained show that the effects of 
viscosity are stabilising and also the increase in the growth rates as — n increases is not 
sufficient to counter the stabilising effects of viscosity which may be quite significant at low 
Reynolds numbers. 

In the past few years it has also been established that additional unstable modes exist 
in which viscosity exerts a destabilising influence. These viscous modes were first discussed 
by Khorrami (1991) and later analysed by Duck Sz Khorrami (1992). Also Khorrami (1992) 
looked at the full viscous incompressible equations by numerical techniques and found that 
at large Reynolds numbers, with n = —1, near the upper neutral point, the first and second 
modes cross and switch. In addition a second peak in the growth rate curves appears due 
to viscous effects. 

In this paper, we take cylindrical polar coordinates, (Lr,9, Lx), with the x-axis lying 
along the axis of the vortex (which is taken to be axisymmetric), and L is some streamwise 
scale. We also take the flow far from the vortex centre to be directed along the x direction. 
The velocity field is written as U * u = U*(u,v,w), the fluid density is p^p, temperature 
T^T, first and second coefficients of viscosity p%aP, p^ A respectively, and pressure p^U* 2 p. 
Here superscript asterisk and subscript oo respectively denote dimensional and freestream 
variables. U* is a velocity scale, defined by 


U 


c oU^ 

Sxul^Rtoc 


log(xi?eoo) + 


8xi /^ 0 2 Re 00 ’ 


( 1 ) 


where is the freestream velocity, = p^/ p^ is the kinematic viscosity, Co is a constant 
and Reoo is the freestream Reynolds number defined by 
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We have a flow Mach number given by 
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where 7 is the ratio of specific heats and R* the gas constant. The freestream Mach number 
is given by 
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U* 

W M ' 


(4) 


where we expect Moo >■ M. Also we define the Prandtl number to be 


a = 


H*c p 


(5) 


where k* is the thermal conductivity of the fluid and c p is the specific heat at constant 
pressure. We further assume a perfect fluid. 

The flow of a swirling wake in the fax downstream limit was first discussed by Batchelor 
(1964) where the flow was taken to be incompressible and a similarity solution was sought. 
This work was extended into the compressible regime by Stott & Duck (1994) who showed 
the following non-dimensional profile was appropriate to leading order 
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where q is a measure of swirl, and 
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It is (6) that we shall be adopting as the basic flow profile for the stability study considered 
in the main part of this paper. We also note that the flow Reynolds number is defined by 


Re = 


PtoU;r, 


( 8 ) 


Consider now the stability equations of the general basic state; taking small perturba- 
tions, we write 


+<*> 


(u,v,w,p, T, p) = (U(r),0,W(r),po,To) 

F(r), iG(r), H(r ), r(r), T(r)l e ^(-ct)/r.+n9) + 
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(9) 


where 8 <C 1 and c = c r + C{ is the complex wavespeed. We assume that p*, the first 
coefficient of viscosity, is solely a function of the temperature T*. In addition we also take 
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the bulk viscosity to be zero. If we now substitute equations (9) and the above assumptions 
into the full Navier-Stokes equations, and consider only terms of 0(6), we obtain 
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Also note that for this study we choose the Prandtl number to be unity and for the purposes 
of the results presented later the ratio of specific heats is taken to be 1.4, although of course 
there is no conceptual difficulty in choosing other values. 

The boundary conditions are given by 


F'(0) = G(0) = H{ 0) = P'(0) = T'(O) = r'(0) = 0 
F(0) = G(0) + H( 0) = P(0) = T(0) = r( 0) = 0 
P(0) = G(0) = H{ 0) = P(0) = r(0 ) = t( 0) = 0 


for 

o' 

II 
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for 

W = i, 

(12) 

for 

Ini > 1, 
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F, G, H,P,T,t -»0 as r — » oo. (13) 

Unlike the inviscid analysis of Stott &: Duck (1994), we are unable to reduce the equations 
further, and thus (in the temporal problem), (10) defines the eigenvalue problem for c. 

In section 2 we develop an asymptotic theory for large wavenumbers when viscous 
effects are included in the analysis. In section 3 it is shown that for sufficiently large Mach 
numbers centre modes are found and consideration is made of the effect of viscosity on these 
modes; indeed, the authors are unaware of previous analysis of the effects of viscosity on 
any kind of centre modes. In section 4 we present and discuss results obtained from the 
analysis of sections 2 and 3. 


2. Analysis for large azimuthal wavenumbers. M=0(|n|*) 

It was shown in the inviscid work of Stott & Duck (1994) that the crucial regime of Mach 
number is when M = 0(|n| * ). In this case there is an important change in the character of 
the inviscid instability. The same is expected to be true in the viscous case, if conditions are 
appropriately chosen. The M = 0(|n|?) regime results also encompass, to a large extent, 
many of the results of smaller orders of M. We then write M = |n|^A/, M = 0(1), 
together with a = no, and the wavespeed is then expected to develop in the form 


c = Co H 1~ 

n 


'Q- 


It is found that the crucial order of the Reynolds number (defined in (8)) is given by 
Re = where /i = 0(1) is the (imaginary) viscosity parameter. Although viscosity is a 

function of temperature, since the temperature is (approximately) constant, then so too is 

A- 

Following the inviscid analysis of Stott & Duck (1994), we fully expect the eigenso- 
lutions to be concentrated close to the point where y>o( r o) = 0, together with <^ 0 ( r o) = 0, 
where 

W 

</>o( r ) = at(U — c 0 ) + — . (15) 

The relevant radial scale is then 

R = (r - r 0 )|n|U (16) 

On the R = 0(1) scale, <p develops in the form 


v = v, + S +0 (l^)- 
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Talcing the 0(1) terms in the governing equations (10) leads to 


r W + « 2 r 0 2 )p 

Ui - ar 0 W o 

(19a) 

Ho = — 

(19b) 

P 

*■ 0 — I ^0, 

U'o - ar 0 Wi 

(19c) 

r 0 /0?(1 + « 2 r 2 ) x 

r ° W'„ ! W i -dr 0 W: + 

(19d) 

where + ^-(l + a 2 ^), and then (p\ is given by the solution of the following quadratic 

equation 
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and hence from this we may deduce ci (since (p \ — — aci). 
of 0(|n|“*), we have that 

Continuing, by equating terms 
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where A\ is an undetermined constant. On account of (20), a solution only exists if the 
following equation is satisfied 
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Equation (22) is of the form of Weber’s equation (after the application of the transformation 

£ = 2? A 4 R). Therefore, in order to satisfy the boundary conditions we must have solutions 
of the form 


_+ 

Ao 

We note that setting /}, the viscosity parameter, to zero reduces the problem to its inviscid 
counterpart as considered by Stott & Duck (1994). Also we see that when \jx\ becomes 
sufficiently large, all disturbances become neutrally stable. 

As in the case of inviscid flow, this regime of M indicates a watershed in our analysis. 
When M is not large, both roots of c\ (obtained from equation (20)) are complex and 
hence the flow is unstable. However, increasing M eventually causes c\ to become purely 
real, which implies C 2 also becomes real. Consequently to this order of approximation, for 
sufficiently large M , all disturbances become neutrally stable. The effects of viscosity on 
this phenomenon and other results obtained from this section will be discussed again in 
section 4. 

We expect, as in the inviscid case, the emergence of centre modes as M becomes large, 
and it is this class of modes that comes under consideration in the next section. 

3. Centre Modes. M=0(|n|) 

We anticipate in this regime, as in the inviscid case, that centre modes (be. modes whose 
disturbances are primarily dictated by conditions close to the axis of the flow) will exist. We 
set M = |n|M, M = 0(1) and assume that disturbances will develop in a similar way to 
those found by Stott & Duck (1994). Thus we expect the flow to be divided into five regions, 
two regions within r = 0(1) (one above, the other below a critical layer), r = 0(|n|~ 5 ) and 
finally r = 0(|n| -1 ). Also, as in the purely inviscid case, we expect the complex wavespeed 
to develop as 

c = Co H — — + . . . , (25) 

n 
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and therefore 


¥> = rvfo H h . . . 

n 

= n{«(f/ - c 0 ) + y } + {-—- + |^o(0) r2 } + ••• • 

For self-consistency, it turns out that we must have <po (r = 0) = 0, t.e. 

c 0 = l + i. 
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(26) 


(27) 


We first consider the inner lengthscale, R = r|n| = 0(1). The other regimes, mentioned 
earlier, will be discussed briefly later, since details are identical to the inviscid case. 

The solution is somewhat complicated by the fact that (as seen in the inviscid solution 
of Stott & Duck 1994) the solution may be regarded as of WKB/multiple scales type. In the 
case of the inviscid solution, it is clear that R derivatives of physical quantities are generally 
0(n). This will then certainly also be the case here. We write 


(P,F,G,H,T,t) =(-Po,F o ,G 0 ,Ho,To,t o )+-(-P u F u G u H 1 ,T u t 1 ) + ... , 
n n n 

The leading order terms in the momentum and continuity equations then give 
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At the next order in n, the following results are obtained: 
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We have implicitly assumed Re = 0(n 5 ); it is this scaling that leads to the key regime 
where viscosity substantially influences the stability of these flows. We therefore write 
A*i = 7 ^. However the problem for H\ =0(1) leads to an eighth-order system, almost of 
comparable computational difficulty to that of the full system as represented by (10). In 
order to make further progress without substantial computational effort, instead, we now 
assume the limit of diminishing viscosity, i.e. fi\ — + 0, (but with |/xi| >> |^|). This enables 
us to write 


(P 0 ,P i ,Go,G l ,H 0 ,H u F 0 ,T 0 ,r 0 ,c 1 ) = (Poo,Pio ,Ooo, Gio,Hoo,Hio, Poo, Too, Too^io) 

+/*i(Poi , P11, 001 , On, Hoi , H\\ , F01, Toi, Toi, cn) + . . . . ( 33 ) 

It turns out, predictably, that Poo is then the leading order inviscid solution, as found 
by Stott &: Duck (1994), and is described by 


£{Poo} = P 00 RR + n2 [ ~~~T ~ - ^-]Poo = 0 . 
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After some algebra, the next order equation (in fi\) is given by 
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Since the solution must be bounded as R — > 0, we have A 40 = A 41 = 0. In the limit 
as R — ► oo, we find 
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Ri(A* 10 R 3 + A* 20 R->), 




Poi ^ RHA* U R° + A* 2 ,R- s ), 


(41) 


Vo (°) 


This solution must be matched to a layer described by f = r|n|? = 0(1), for which 
the following solution can be written 


= r* (<?,/„( 


,2 » / f 2 |^| ;; 2 „ 2 ~ 2 


q 2 M |n|r 2 
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) exp (i / 
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where J„ and K„ are modified Bessel functions with v = Note that v (and also s) is 
imaginary. A relationship between C\ and Ci is determined from a consideration of the 
r = 0(1) region, which is inviscid in nature, and leads to a WKB-type solution, giving the 
result 


C\ — —2ittC2 exp 


l/Xi- 


-j-»2 2 W A W 
-M y> 2 + 


4r 2 
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where r*o is determined from 
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Matching asymptotically, and taking the limit |n| — > oo, leads to the conclusion (on account 
of (43)) 


-^io — -4-n — ■ • • — 0. 


(45) 
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Equation (34) may be solved using hypergeometric functions (Stott & Duck 1994), and this 
leads to 


2n 2 q{\ + qoc) 

WWW * 1 + 1 J ’ 


(46) 


where ft = |[|n| + s] + N , N = 0, 1, 2, 3 . . . Alternatively a WKB solution would also appear 
to be in order here, although this is not straightforward due to the multitude of critical 
points. Since the system (35) depends on a largely numerical approach, we chose to adopt a 
numerical approach for the problem overall, including (34), treating n as a finite (but large) 
numerical parameter, using a fourth order Runge-Kutta scheme. 

This analysis differs somewhat from that of Stewartson & Brown (1985), who looked 
at near neutral centre modes. In their work s was real (and so too, consequently was Cio), 
and it was necessary to seek higher order terms in order to obtain the imaginary part of cj . 
In the following section, we consider some numerical results arising from the analysis of this 
and the previous section. 


4. Conclusions 

In the work described in the previous two sections we find that as in the corresponding 
inviscid analysis, the Mach number proves to be crucially dependent on the order of the 
wavenumber and we are able to develop the analysis in a similar way to the previous inviscid 
study. Note that for all results referred to in this section the swirl parameter q is taken to 
be 0.8. 

Khorrami (1992) noted that modes evaluated above Reynolds number as low as 10000, 
attain their maximum growth rates and beyond this point appear to be largely independent 
of the Reynolds number. Figure 1 shows a plot of growth rate against Reynolds number 
for n = —15 with a = 7.5 with M = |n |2 (for the first mode i.e. that corresponding to 
m = 0 in (24)). We observe from these results that stabilisation occurs for sufficiently low 
Reynolds number; this is a generic feature of this problem (and indeed also occurs in the 
incompressible case, generally). 

Figure 2 shows the neutral curve (ct against Reynolds number) for the case n = — 15 
with M = 0(|n|s) (again for the first mode), and confirms the remarks regarding the effects 
of stabilisation above. 

Figure 3 shows the variation of the growth rates with Mach number for the most 
unstable mode for varying Reynolds numbers and n = —15, with M = 0(|n| 2 ). In this 
figure it can be seen that for all Reynolds numbers considered the effect of increasing Mach 
number is generally stabilising and for sufficiently large M modes are completely stabilised. 
Also as the Reynolds number is increased the value of the Mach number at which the modes 
become stable decreases. 

We now turn to consider the effects of viscosity on centre modes. Figures 4a-4c show 
growth rate variations for the case M = 15, n = —15, q = 0.8, <7 = 1, for Re — 00 (fig. 
4a), Re = 2.5 x 10 6 (fig. 4b) and Re = 10 6 (fig. 4c). These figures show the first three 
modes, with the first mode ( N = 0) shown as a solidus, the second mode ( N = 1) shown 
as the larger dashed line, and the third mode ( N = 2) shown as the smaller dashed line. 
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The first of these figures, namely fig. 4a corresponds with the inviscid limit, as considered 
by Stott & Duck (1994), and agrees very favourably with these previous results, obtained 
by (46). These results also show the trend observed by Stott & Duck (1994), namely that 
near the lower neutral point (a — > 2 in this case), that the higher modes do have slightly 
larger growth rates than corresponding lower modes. Figures 4b, 4c do generally indicate 
that viscosity plays a stabilising effect on all three modes considered. Indeed, these observed 
trends are, again, quite generic, and are observed in general calculations of this problem. 
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Figure 1: Variation of growth rates against Reynolds number for M — |u| ^ , n — —15, 
a = 7.5. 
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effects become important at much larger Reynolds numbers than for the first class of disturbance. 
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